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A SIMPLE  MODEL  OF  A LINUS  FUSION  SYSTEM 


WITH  A THICK,  COMPRESSIBLE,  RESISTIVE  LINER 


I.  Introduction 

A simple  model  of  the  dynamic  behavior  of  a cylindrical  Incompres- 
sible, perfectly-conducting  liner  Imploding  onto  a g = 1 plasma  was  de- 
veloped early  In  the  LINUS  program^  In  order  to  derive  the  basic  scaling 
laws  for  power-producing  systems.  The  principal  relationship  obtained 


Q = Bt,roP%(b,  To) 


where  Q = nuclear  energy  released  In  one  cycle/plasma  energy  at  peak 
compression; 

Bg  = buffer  magnetic  field  between  liner  and  plasma  at  peak  com- 
pression (asstimed  of  negligible  extent); 

Tq  = minimum  liner  radixis; 
p = density  of  liner  material; 

Tq  = plasma  temperatvire  at  peak  compression; 
b = a dimensionless  thickness  parameter,  defined  by  A = irb^roj 
where  A Is  the  cross-sectional  area  of  the  liner. 

The  function  P(b,  Tg)  was  evaluated  by  integrating  the  thermonu- 
clesir  reaction  rate  over  a compression-expansion  cycle  (assumed  symmet- 
rical). It  has  a shallow  maximum  at  Tg  = 20  Key,  and  varies  with  b as 

shown  in  Figure  1. 

Note:  Manuscript  submitted  March  24,  1977. 
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An  iB^Jortant  featijre  of  the  model  was  the  description  of  the  tra- 
jectory in  time  units  of  r /v  , where 

O 00 


is  the  velocity  of  the  liner  if  allowed  to  expand  to  infinity.  The 
"dwell  time"  around  peak  compression  during  which  reactions  occur  is 
then  given  by  ^ f(b). 

09 

In  Reference  1,  the  effects  of  liner  resistivity  and  compressibil- 
ity were  briefly  considered.  The  effect  of  resistivity  is  to  allow  the 
buffer  magnetic  field  to  penetrate  into  the  liner  and  thus  seme  of  the 
liner  energy  is  dissipated  that  otherwise  would  go  into  plasma  compres- 
sion. Provided  that  the  skin  depth  is  small,  the  effect  of  resistivity 
may  be  treated  as  a perturbation,  and  the  treatment  of  Reference  1 is 
still  basically  valid,  although  more  detailed  calculations  in  which  the 
diffusion  is  incorporated  self-consistently  into  the  liner  dynamics 
would  lead,  to  greater  accuracy. 

The  effect  of  compressibility  was  also  treated  as  a perturbation  by 

asstming  that  all  elements  of  the  liner  come  to  rest  at  peak  compression 

and  calculating  the  ratio  E /E  of  the  energy  stored  in  compression  of 

c o 

liner  material  to  the  energy  in  the  compressed  plasma.  A characteristic 
magnetic  field  was  defined  by  *=  SttPc^®,  where  c^  is  the  speed  of 

sound  in  the  liner  material.  Table  1 sho^-ra  values  of  B^  for  some  can- 
didate liner  materials.  The  fraction  of  energy  in  material  compression 
was  then  given  by 
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The  function  g is  shown  in  Figure  2.  For  the  thin  liners  (b  ~ 3)  that 
were  being  considered  in  the  early  stages  of  the  LINUS  program,  the  as- 
sumptions of  the  above  treatment  were  justified  provided  Bq  < Bj . 

The  situation  has  changed  drastically  now  that  we  are  considering 
captive  liquid  liners  which  will  be  very  thick  (b  ~ 30-4o)  at  peak  com- 
pression. Initial  results  from  a one-dimensional  fluid  code®  show  that 
the  inside  surface  comes  to  rest  while  the  outside  surface  is  still  mov- 
ing; the  primary  effect  of  compressibility  is  thus  to  modify  the  dynam- 
ics due  to  the  finite  speed  of  so\md  and  to  inhibit  the  transfer  of  en- 
ergy from  the  outer  regions  of  the  liner  to  the  inner  regions  on  the 
time  scale  of  the  turnaround.  The  res\iit  is  that  the  inside  surface  of 
the  liner  turns  roiond  at  a larger  radius  and  produces  a lower  peak  plas- 
ma pressure  than  would  have  been  anticipated  by  the  incompressible  model. 

The  purpose  of  this  note  is  to  describe  a simple  model  which  en- 
ables the  actual,  turnaround  radivis  r^  to  be  derived  in  terms  of  the 
tTimaround  radius  r^  anticipated  by  the  incompressible  model.  The  pres- 
ent model  assumes  that  initieilly  the  liner  motion  can  be  predicted  by 

the  incompressible  dynamics  used  previously,  but  that  at  some  point  in 

I 
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1 the  compression  cycle  the  liner  can  be  divided  into  two  regions,  com- 

I 

pletely  decoupled  from  each  other.  The  outer  region  has  no  further  ef- 
fect on  the  inner  region  which  executes  a symmetrical  compression- 
expansion  trajectory,  again  determined  by  incompressible  dynamics.  The 
determination  of  the  point  of  separation  involves  an  adjustable  para- 
. meter  which  may  be  used  to  fit  the  model  to  more  exact  numerical  calcu- 

) 

; lations.  However,  it  will  be  shr'wn  that  the  results  are  relatively  in- 

[ sensitive  to  this  parameter,  which  leads  us  to  believe  that  the  model  is 
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sufficiently  accurate  to  be  used  in  scaling  studies  and  conceptual  de- 
signs of  LINUS  systems. 

II.  The  Model 

The  central  feature  of  the  model  is  the  determination  of  the  point 
of  separation  of  the  inner  and  outer  regions  of  the  liner.  The  liner 
position  at  the  point  of  separation  is  shown  in  Figure  3*  The  radius  of 
the  inside  surface  is  the  radius  of  the  outside  surface  is  and 
the  radius  of  the  sxirface  of  separation  is  r.  The  inner  shell  of  the 
liner  will  eventually  come  to  rest  at  radius  r^,  which  is  of  course 
greater  than  the  raditis  rj,  at  which  the  entire  liner  would  have  come  to 
rest  on  the  incompressible  model. 

Using  the  dimensionless  thickness  parameter  b,  defined  in  the  In- 
troduction, we  have 

If  = a^r%  + bsps  (4) 


Similarly,  we  may  define  a thickness  parameter  b^  for  the  inner 
shell  only 


? = 


(5) 


The  energy  that  will  be  delivered  to  the  plasma  when  the  inner 
shell  comes  to  rest  is 


where  Eg  is  the  total  energy  initially  in  the  liner 
Eg  = irrr|bapv2  . 
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(6) 

(7) 


hill  ill  lift 


Continuing  the  convention  that  starred  quantities  represent  the 


values  actually  achieved  by  the  compressible  fluid  we  find 


> — =( — T 

Po  ~Y^)  ’ Bo  \r^/ 


The  energy  that  has  been  delivered  to  the  plasma  at  the  point  of  separa- 


tion is 


2(V-1) 


The  energy  remaining  to  be  delivered  is 


The  energy  in  the  inner  shell  of  the  liner  is 


E,  . (E„  -E,) 

I / , \2(V-1)\  / 


Equating  and  Eg , and  putting  S = , we  obtain 

^0 


1)  (l  . = ((aS)=<''-y.  1)  in  (l  . 

This  equation  may  be  used  to  obtain  b^  as  a function  of  S,  with  a and  b 
as  parameters. 

We  now  define  an  energy  delivery  time  as  the  time  taien  for  the  in- 
ner surface  to  go  from  ar^  to  r^ 

tj  = ^ F(b^,a,  V)  (13) 

* 
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where  F is  obtained  from  the  rather  simple  calculation  of  the  incom- 
pressible trajectories  described  in  Ref.  (1).  Since  only  the  inner 
shell  is  involved,  we  have  from  (2) 


I 


V 2 
* 


_2p*_ 

(Y-l)b|p 


(14) 


We  also  define  an  energy  communication  time,  as  the  time  for  a 
sound  wave  to  travel  from  the  outside  to  the  inside  of  the  inner  shell: 

ta  = (r  - ((1  + - 1)/Cg  (15) 

The  assumption  of  incompressible  motion  of  the  inner  shell  is  that 
tg  i tj  , We  now  choose  to  define  the  point  of  separation  by  taking 
tg  = tj^ , Introducing  a new  parameter 


C = Po/Pc|  = B|/B? 


(16) 


which,  as  described  in  the  Introduction,  may  be  used  as  a measure  of  the 
degree  to  which  compressibility  must  be  considered,  we  obtain  from  equa- 
tions (8)  and  (13)  through  (16) 

(V-l)b2  s2ViC(b^,a,  y)  = 2C((1  + b»  )i  - 1)2  (17) 

This  equation  may  be  used  to  obtain  S as  a function  of  b^  , with  a and  Q 
as  parameters.  Given  b,  Q and  a,  we  may  now  use  equations  (12)  and  (17) 
to  determine  S and  b^  by  graphical  methods,  b and  ^ characterize  the 
aimed- for  conditions  in  terms  of  the  earlier  incompressible  model.  S and 
b^  characterize  the  achieved  conditions  predicted  by  the  current  model, 
and  all  other  relevant  quantities  may  be  derived  from  them. 
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Ill,  Choice  of  a 
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The  point  of  separation  may  be  defined  in  terms  of  a fraction 

^ _ energy  remaining  to  be  delivered  at  point  of  separation 
~ total  energy  delivered  at  peak  compression 


which  is  related  to  the  parameter  a by 


f = 1 


a 


-a(Y-l) 


(18) 


As  previously  stated,  a is  an  adjustable  parameter  which  may  be  used  to 
fit  the  results  of  the  present  model  to  the  results  of  more  exact  hydro- 
dynamic  computations.  A particular  case  for  which  computed  results  are 
available^  is  b = 4o,  g = 0.25.  We  take  these  initial  conditions  and 
use  the  present  model  to  derive  S and  b^  using  four  values  of  a,  namely 
l.l8,  1.35»  1.68,  and  2.28,  which  correspond  to  values  of  f of  0.20, 
0.33»  0-50,  and  0.6?  respectively  when  Y is  taken  as  5/3*  The  results 
are  shown  in  Table  2.  Note  that  for  f < 0.5»  S varies  by  less  than  1'^ 
from  the  value  1.32.  b*  increases  steadily  as  f increases,  as  might 
be  anticipated  from  simple  physical  considerations.  The  computed  value 
of  S is  1.57,  which  cannot  be  reconciled  to  our  model  simply  by  adjust- 
ing a.  As  will  be  shown  in  the  next  section,  this  discrepancy  is  re- 
solved when  the  energy  stored  in  compression  of  the  liner  material  is 
taken  into  account . 

It  is  difficult  to  estimate  b*  from  the  computed  trajectory,  since 
the  latter  is  unsymmetrical  with  respect  to  the  point  of  turnaround; 
however,  b»-  ~ 10  provides  a reasonably  good  fit.  In  all  subsequent 
calculations,  we  shall  therefore  give  0!  the  value  of  1.68,  thus  defining 
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the  point  of  separation  as  the  point  where  50^  of  the  final  energy  has 
been  delivered  to  the  plasma.  Table  3 summarizes  the  results  of  20 
cases  which  cover  the  parameter  range  of  practical  interest  and  ift 
addition  to  S and  b^,  lists  some  other  derived  quantities. 

IV.  Correction  for  Material  Compression 

In  the  preceding  discussion,  the  effect  of  compressibility  was 
manifest  only  through  the  finite  speed  of  sound  which  brought  about 
the  decoupling  of  the  inner  and  outer  shells  of  the  liner.  Since  the 
inner  shell  is  itself  quite  thick  (4  < b^  < 12)  it  is  necessary  to  take 
into  account  the  energy  stored  in  the  compression  of  the  liner  material 
when  calculating  the  motion  of  the  inner  shell.  From  equation  (3)  we 
find  that  the  fraction  of  the  inner  shell  energy  stored  in  material 
compression  when  the  shell  comes  to  rest  is 

= S(l>.)CP»/Po  (19) 

where  g is  the  function  shown  in  Fig.  2.  is  tabulated  on  the 

bottom  line  of  Table  3;  it  amounts  to  ~ 20-25^  and  thus  repre- 
sents a small  but  significant  perturbation.  Although,  strictly  speak- 
ing, the  energy  of  compression  should  be  incorporated  self-consistently 
in  the  calculation  of  the  motion,  it  is  in  keeping  with  the  simple 
nature  of  our  model  to  consider  the  effect  only  as  a correction  to  the 
calculations  already  made.  Subtracting  the  energy  of  compression  from 
the  energy  of  the  plasma  results  in  an  increase  in  S by  a factor 
(1  - E^/E^)  with  a corresponding  effect  on  the  derived  quantities. 
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The  corrected  results  are  displayed  in  Table  4.  Note  that  for  our 
test  case,  b = 40,  Q = 0.25,  the  model  now  predicts  S = 1.55»  compared 
to  the  computed  value  of  1.57.  This  remarkably  good  agreement,  although 
probably  fortuitous,  gives  one  some  confidence  in  the  validity  of  the 
present  model. 


V.  Thermonuclear  Reactions 

The  thermonuclear  yield  is  obtained  by  applying  the  analysis  of 
Ref.  1 to  a plasma  compressed  by  the  inner  shell  only.  If  Q is  the 
ratio  of  the  fusion  energy  yield  to  the  energy  initially  in  the 
entire  liner,  and  we  neglect  the  energy  in  compression  of  liner  material 
we  have 


Q » 


T 

o 


(20) 


where  F is  the  function  shown  in  Fig.  1,  Bi  is  the  characteristic  field 

for  the  liner  material  given  in  Table  1,  and  the  primes  are  used  to 

indicate  uncorrected  values  of  the  quantities. 

The  correction  for  material  compression  is  made  by  taking  the 

dwell  time  as  unchanged,  but  multiplying  the  radius  by  (1  - E^/E^) 

1-5 

and  the  plasma  density  by  (1  - E^/E^)  . Thus 


« = (f)  (r ) c®'  O’,  1 V* 


We  now  take  T^  = 20  keV  and  construct  a new  function  from  equation  21, 


= Q/Bir^p^  (22) 

is  also  tabulated  in  Table  4.  The  units  are  the  same  as  in  Fig.  1. 
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Given  Q,  b and  the  material  constant  (Table  1),  enables  the 

compressed  plasma  radius  r^  to  be  determined  for  a given  Q,  and  hence, 
through  the  other  quantities  listed  in  Table  the  scale  of  the  whole 
system. 

An  important  quantity  is  the  total  (initial)  liner  energy  per  unit 
length  given  by 

2 

For  our  present  purpose  it  is  convenient  to  express  the  energy  in  terms 
of  a new  quantity: 


E(g,b)  = E^p/Q^  (p  in  gm.cm  E^  in  MJ.m"^) 


(2h) 


values  of  which  are  included  in  Table 

As  will  be  apparent  in  the  next  section,  the  time-scale  for  energy 
delivery  t^  is  an  important  quantity,  defined  by  equations  (IJ)  and 
(lij-) . With  a = 1.682  and  y m 5/5  we  have 


(25) 


where,  as  previously,  we  assume  that  this  time  is  not  altered  when 
material  coaq>ression  is  taken  into  account.  We  may  now  construct  a 
dimensionless  time,  as  follows 


t - ticg/r^ 


(26) 


values  of  which  are  given  in  the  last  part  of  Table  4. 

A number  of  Interesting  conclusions  can  be  drawn  from  the  results 
shown  in  Table  4.  We  note  first  that  over  the  entire  range  of  para- 
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meters  studied  (b  = 10-50,  ^ a 0.1-2.5)»  varies  by  less  than  ±25^. 
This  means  that,  given  the  Q of  the  system  and  the  liner  material,  the 
compressed  plasma  radius  r^  is  almost  Independent  of  Q and  b.  This  is 
in  direct  contrast  to  the  prediction  of  the  incompressible  model, 
namely  that  r^  <*  which  previously  led  us  to  believe  that  the 

highest  magnetic  fields  were  desirable  to  reduce  the  radial  dimension 
of  a LINUS  fusion  system.  The  incompressible  model  also  predicted 

A 

that  the  energy  per  unit  length  E was  independent  of  B , but  the  com- 

o 

A 

pressible  model  indicates  that  E increases  as  Q increases. 

Clearly  as  ^ is  reduced  there  must  come  a point  Q = Qib)  at  which 
the  incompressible  model  becomes  applicable;  this  point  can  be  found 
by  putting  S * 1,  b * b and  £ » C in  equation  (17)  and  solving  for  Q. 

A ^ 

Results  are  given  in  Table  5»  which  also  lists  and  E at  C 
Although  at  this  point  soaid  transit  effects  are  not  important,  material 
compression  must  still  be  considered  since  although  C is  small,  b s b^ 
is  large  and  E^/E^  ~ 10-155^.  This  correction  has  been  applied  in 

A 

calculating  P^  and  E.  The  principal  conclusion  that  can  be  drawn  from 
Table  5 is  that  in  most  cases  of  practical  interest  Q >%  and  so  the 
compressible  model  must  always  be  used. 

We  conclude  this  section  with  a numerical  example.  Consider  a 
'breakeven'  system  (Q  = 1)  with  a lead  liner,  and  take  ^ a 0.25,  b = 40 
as  operating  parameters.  From  Table  4,  we  find  r^  • 3 cm  and  E^  a 
23  MJjn  The  peak  magnetic  field  = 0.8  MG.  A system  designed  to 
achieve  fusion  energy  = plasma  energy  needs  only  Q = E^/E^  a 0.57, 
in  which  case  r^  = 1.7  cm  and  E^  a 7-3  MJ.m  These  figures  should 
be  compared  with  r^  a 1 cm  and  Eq  a l4  MJ.m  ^ derived  in  Ref.  1 for  a 
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Q ss  1 system  with  a thin  copper  liner  at  a 3 MG. 


VI.  Effect  of  Finite  Resistivity 

The  effect  of  compressibility  is  to  reduce  the  efficiency  of 
energy  transfer  from  the  liner  to  the  plasma;  the  energy  not  delivered 
is  stored  in  fluid  motion  and  material  compression  and  can,  in  principle, 
be  recovered  mechanically.  On  the  other  hand,  the  finite  resistivity 
of  the  liner  results  in  an  irreversible  loss  of  energy  due  to  resis- 
tive dissipation  of  the  flux-conserving  current  on  the  inside  surface 
of  the  liner,  and  this  energy  appears  as  heat.  As  in  Ref.  1,  we  shall 
treat  this  effect  as  a perturbation,  that  is,  we  shall  assume  6/r^«  1, 
where  6 is  the  skin  depth.  We  assume  that  conditions  at  peak  com- 
pression are  not  changed,  and  that  the  energy  required  to  compensate 
for  the  dissipation  is  obtained  by  increasing  the  initial  energy  of 
the  liner. 

We  are  mainly  interested  in  the  total  energy  lost  in  a complete 
compression-expansion  cycle,  and  so  we  assume  that  the  inside  surface 
of  the  liner  is  subjected  to  a magnetic  pulse 

B = sin  (2n/T)t  (27) 

and  use  the  result  of  Kncepfel^  that  the  energy  skin  depth  at  t = T/2 
is  given  by 

6 = 1.27  (T/1ht^)^  (28) 

where  a is  the  electrical  conductivity  of  the  liner  material  (6  in  cm, 
a in  e .m.u. ) . 


f 


I 


i 
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The  magnetic  pulse  at  the  liner  surface  is  not  in  fact  sinusoidal, 
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but  the  significant  part  of  the  pulse,  which  occurs  between  the  point 
of  separation  and  peak  compression,  can  be  approximated  as  a sine  wave 
with  fair  accuracy,  and  this  approximation  is  in  keeping  with  the 
simple  nature  of  our  model . 

At  the  point  of  separation 

B = B » B ' (1.682)‘®'^3  = 0.420  B^(l  - f (29) 

s *■  * c 

where  f « E /E'.  The  matching  to  the  sine  wave  is  obtained  by  putting 
c c ^ 

ti  * kiT/2rr,  where  ti  is  given  by  equation  (25)  and  k^  is  given  by 

ki  a cos"^(B^/B„)  s cos‘^(0.420  (1  - f (30) 

S ” c 


and  substituting  r^  = Q/P^B^p^  and  Bi  a (8npc®)^10  ® (Bx  in  MG)  we 


obtain 


/plo’®  /i.29Pj\^ 

f.-l— ) 


The  dimensionless  quantity  L is  tabulated  in  Table6.  it  varies  by  less 
than  a factor  of  2 over  the  range  of  parameters  studied. 

As  an  example,  consider  a system  defined  by  Q a 1,  Q m 0.25, 
b a 4o  and  take  a lead  liner  (p  « 11 -5  gm. cm  ®)  with  a thin  lithium 
layer  (a  =2.10  ®e.m.u.)  on  the  inside  to  carry  the  skin  current. Then 
6/r^j,  =»  0.2,  and  so  the  effect  of  resistivity  will  be  small,  but  signifi- 
cant. The  energy  lost  is  E (2/3) (25/r„) ; the  factor  2/3  arises  because 
the  lost  energy  is  in  magnetic  field  rather  than  plasma. 


1 


VII.  LINUS  Reactors 


We  make  optimistic  assumptions  regarding  the  plasma,  namely  that  it 
is  P ■ 1,  perfectly  contained  and  that  C = 0.2  (all  a-particles  trap- 
ped) . We  take  a composite  liner  with  the  density  of  lead  and  the  con- 
ductivity of  lithium  and,  with  b a 40,  we  solve  equation  (3^)  for  Q as 
a function  of  ^ and  the  product  TIiTla*  The  results  are  shown  in  Table  7. 
We  see  that  for  systems  with  very  efficient  fluid  transfer  (T)iTli  ^ 0.9) , 
a self-sustaining  reactor  cycle  is  possible  with  Q < 1.  We  note  also 
that  Q decreases  as  Q increases,  but  by  less  than  a factor  of  2 as  ^ 
goes  from  0.1  to  2.5.  As  an  example,  we  choose  ^ = 0.25,  T]iTla=r0.9, 
giving  Q a 1.25  and  calculate  the  principal  parameters  of  a LINUS  re- 
actor. The  results  are  shown  in  Table  6.  This  choice  of  parameters 
conveniently  leads  to  a liner  thickness  of  ~ 100  cm  at  peak  compression, 
which  is  a desirable  thickness  for  neutron  absorption  and  tritium 
breeding. 

As  we  have  seen  throughout  this  report,  the  system  parameters 
do  not  in  general  depend  strongly  on  the  choice  of  An  exception  is 
the  energy  per  unit  length  E^,  and  through  this,  the  pressure  required 
to  drive  the  liner.  Consider  a captive  liquid  liner^>^  driven  by  a 
constant  external  pressure  Since  with  this  arrangement,  the  dis- 

placement of  the  liner  equals  its  volume,  we  may  equate  the  work  done 
by  Pj^  to  the  liner  energy  and  obtain 


where  R is  the  additional  fraction  of  energy  required  for  rotational 
stabilization.^  Taking  R = 1,  Bj.  = 3-2  MG  (lead),  y = 5/5  we  have 

= 1.22C/b2  Mbars  (36) 

with  b « 40  and  Q = 0.25,  equation  (36)  gives  =*  I90  bars  *»  2800  psi, 
which  is  somewhat  less  than  the  steam  pressures  (m  3300  psi)  used  in 
advanced  thermal  power  plants.®  Conversely,  it  might  be  argued  that 
the  engineering  limit  on  driving  pressure  restricts  the  choice  of  Q to 

< 0.3. 

Finally  it  should  be  said  that  the  principal  shortcoming  of  the 
model  presented  here  is  the  assumption  of  a perfectly  confined  ^ s 1 
plasma.  To  restrict  end  loss  in  a system  of  reasonable  length  it 
appears  that  a containment  geometry  with  closed  magnetic  field  lines 
will  be  necessary^  and  such  systems  generally  have  p ^ 0.5.  A dis- 
cussion of  these  systems  is  beyond  the  scope  of  this  report  but  the 
effect  on  our  results  can  at  least  be  estimated. 

The  effect  of  p < 1 is  to  require  that  P (b,  T^) , equation  (1)  and 
Fig.  1,  should  be  multiplied  by  a factor  s,  where  s ~ < 1.  Then 

from  equations  (20)  to  (22)  we  see  that  must  be  multiplied  by  the 
same  factor.  Since  P^  occurs  again  in  equation  (32),  A must  be 
multiplied  by  s^,  which  in  turn  affects  the  value  of  Q derived  from 
equation  (34).  If  11  ilia  « 1,  Q <*  and  hence  Q must  be  multiplied  by 

s^^®.  Since  r^  * Q/P^,  r^  must  be  multiplied  by 

The  principal  consequence  of  p < 1 is  therefore  to  increase  the 
size  of  the  system;  some  preliminary  calculations  by  the  author  give 
s ~ 0.25  for  a typical  closed-field  configuration.  This  would  increase 
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by  2.5  and  the  reactor  outlined  in  Table  8 would  become  rather 
cumbersome.  This  study  therefore  indicates  the  need  for  plasma  con- 
finement with  ^ as  close  as  possible  to  unity. 
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TABLE  3:  Results  uncorrected  for  material  compression 


TABLE  4 ; Results  corrected  for  material  compression 


(a)  Table  of  S = r^/r^ 


10 

20 

50 

4o 

50 

0.1 

1.21 

1.51 

1.55 

1.58 

l.4o 

0.25 

1.55 

1.44 

1.50 

1.55 

1.54 

1.0 

1.64 

1.76 

1.81 

1.87 

1.89 

2.5 

1.90 

2.05 

2.11 

2.15 

2.18 

<b) 

Table  of 

C\ 

10 

20 

50 

4o 

50 

0.1 

8.4 

10.1 

10.7 

11.5 

11.7 

0.25 

5.8 

6.8 

7.4 

7.8 

8.1 

1.0 

5.5 

4.1 

4.4 

4.7 

4.8 

2.5 

2.6 

5.1 

5.5 

5.5 

5.6 

(c)  Table  of  B^/B^ 


10 


20 


50 


40 


50 


1 


TABLE  (Con’t.) 


(d) 

Table  of 

P*/Po 

I 

10 

20 

30 

4o 

50 

0.1 

0.53 

o.4i 

0.37 

0.54 

0.33 

0.25 

0.37 

0.30 

0.26 

0.24 

0.24 

1.0 

0.19 

0.15 

0.l4 

0.12 

0.12 

2.5 

0.12 

0.09i 

0.083 

0.080 

•j 

0.074 

■1 

(e) 

Table  of 

* 0 

1 

1 

10 

20 

30 

4o 

i 

30 

O.J 

0.78 

0.70 

0.67 

0.65 

0.64  j 

0.25 

0.67 

0.61 

0.58 

0.57 

0.56  1 

1.0 

0.52 

0.47 

0.45 

0.45 

0.45  ^ 

TABLE  4 (Con’t.) 


(g)  Table  of 


X 

10 

20 

30 

0.1 

0.037 

0.030 

0.029 

0.25 

0.038 

0.033 

0.030 

1.0 

0.036 

0.051 

0.029 

2.5 

0.032 

0.028 

0.026 

(h)  Table 

of  E 

10 

20 

30 

0.1 

9.1  10^ 

1.2  10^ 

1.2  10^ 

0.25 

1.8  10^ 

2.1  10^ 

2.5  10^ 

1.0 

5.4  10^ 

6.3  10^ 

6.8  102 

1.4  10^  1.6  10' 

(i)  Table  of  t 


40 


0.027 

0.028 

0.027 

0.025 


50 


0.026 

0.028 

0.026 

0.024 


I'l 

il 


4o 


1.5  10^ 

2.5  10® 
7.4  10^ 
1.6  10^ 


50 


1.4  10® 

2.5  10^ 
7.8  10^ 
1.7  10® 


'i 


; ' 

f 

h 

I* 


TABLE  5: 


Limits  of  incompressible  model 


b 

10 

20 

50 

4o 

50 

C 

0.066 

0.020 

0.010 

0.006 

0.004 

0.056 

0.025 

0.020 

0.017 

0.014 

A 

£ 

71 

48 

58 

55 

52 

TABLE  8 


Parameters  of  a LINUS  reactor  operating  with  a 
self-sustained  mechanical  cycle 

Dimensionless  parameters:  Q — 0.25 

b = i+0 
Q = 1.25 
lliTla  =0.9 

Liner  material:  Bi  = 3*2  MG  (lead) 

p = 11.3  g®  cm  ® (lead) 
a = 2.10  ® e.ra. u.  (lithium) 

Derived  quantities:  r^  = ^.15  cm  (compressed  plasma  radius) 

ra  =!  108cm  (liner  outer  radius) 

= 800  kG  (peak  magnetic  field) 

E = 36  MJ.m  ^ (liner  energy,  excluding 

rotation) 

E^  = 45  MJ.m  ^ (fusion  energy) 
ti  = 40  ^sec  (~  ^ dwell  time) 

p^^  2800  psi  (driving  pressure,  including 

rotation) 
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